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Abstract. With few exceptions, known explicit solutions of the curve shortening flow 
(CSE) of a plane curve, can be constructed by classical Lie point symmetry reductions or by 
functional separation of variables. One of the functionally separated solutions is the exact 
curve shortening flow of a closed, convex "oval" -shaped curve and another is the smoothing of 
an initial periodic curve that is close to a square wave. The types of anisotropic evaporation 
coefficient are found for which the evaporation-condensation evolution docs or does not 
have solutions that are analogous to the basic solutions of the CSE, namely the grim reaper 
travelling wave, the homothetic shrinking closed curve and the homothetically expanding 
grain boundary groove. Using equivalence classes of anisotropic diffusion equations, it is 
shown that physical models of evaporation-condensation must have a diffusivity function 
\ that decreases as the inverse square of large slope. Some exact separated solutions are 

{SJ ■ constructed for physically consistent anisotropic diffusion equations. 
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1 Introduction 

^ \ The standard curve shortening flow is a nonlinear evolution by curvature, 

— (x,t) = K(x,t)n(x,t), j(x,0) = 70 (x), (1) 

where k(x, t) is the Euclidean curvature of 7(2;, t) and n(x, t) is the outward unit normal to the 
curve at each point. For the curve-shortening flow parameterized as x ^ (x,y(x,t)), 

c?7 11 
n-e y y t = n-—{x,t), y t 



which describes the curve shortening flow in Cartesian coordinates. This equation arises in the 
practical context of metal surface evolution [24]. More recently, it has been used extensively as 
an isotropic curve-smoothing mechanism in image processing |23[ [25] ; see also [9] . Equation ([1]) 
is well known and has been the subject of several important studies, for instance |16[ 117]. Few 
explicit and exact curve shortening flows of plane curves are known in the literature. Some 
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similarity solutions with prescribed-slope boundary conditions have been constructed paramet- 
rically in terms of integrals of algebraic functions of elementary functions [SJ [2DJ . King [2T] writes 
down a number non-invariant solutions. In addition, there are known to be various self-similar 
rotating spiral and flower- head solutions [H [18] . 

We demonstrate in this paper how for most solutions known to us, Lie symmetries of various 
types or separable coordinate systems play a vital role and we conjecture that all currently 
known solutions and many new solutions can be obtained via a series of higher order constraints 
in the method of functional separation. In Section [2] we briefly review the solutions that can 
be obtained by symmetry reduction. In Section [3] we show how the method of functional 
separation of variables recovers additional interesting explicit, non-self-similar solutions. Some of 
the properties and materials science applications of these solutions are developed in more detail. 
In Section [H an anisotropic version of the second-order Mullins equation is derived for materials 
that include a realistic dependence of evaporation coefficient on surface orientation. Equivalence 
classes for anisotropic evaporation coefficients are constructed under the Euclidean group. Using 
the equivalence classes, in Section[5]we investigate under which types of anisotropy, analogs of the 
standard solutions for isotropic diffusion, do or do not exist. Some examples of exact solutions 
are produced by direct construction and in Section [6] by functional separation of variables. 

2 Symmetry reductions of the curve shortening equation 

One of the most widely used techniques for constructing explicit, exact solutions of nonlinear 
partial differential equations is Lie symmetry reduction. This has been extensively investigated 
in relation to the partial differential equation ([2]) or its derivative forms 

u t = d x [D(u)u x ], D{u) = l/(l + u 2 ) , (3) 

Kt = K 2 K d9 + K 3 , (4) 

R t = de[R- 2 R ]-R~ 1 , (5) 

where u = y x , k is curvature, R = 1/k and 6 is the orientation angle along a convex curve [2]. 
Although the standard nonlinear diffusion equation ([3]) is rarely used in the context of curve 
shortening, it has non-trivial Lie potential symmetries that enable one to construct exact simi- 
larity solutions [3]. The reaction-diffusion equation in standard form ([5]), another equation that 
is rarely used in the context of curve shortening, has been fully classified by Lie point symmetry 
reductions [15]. 

In common with all autonomous nonlinear diffusion equations of second order, the curve 
shortening equation ([2]) is invariant under translations in x, y and t, plus the Boltzmann scaling 
group generated by xd/dx + ydy + 2td/dt: 

(x,y,t) = (e € x,e e y,e 2e t). 

This allows the possibility of two types of scale-invariant similarity solution: the expanding 
solution of the form y/y/t — to = G{x/y/t — to), (t > to) an d the contracting solution of the form 
y/y/to — t = G(x/ '\/to — t), (t < to). Being Euclidean-invariant, ([2]) has the rotation group as an 
additional symmetry. Some interesting exact solutions to this equation, and consequently to ([3j) 
can be constructed by consecutive symmetry reductions [4]; see also [TT] . 

The expanding homothetic similarity solution with initial-boundary conditions y x (0,t) = m, 
y — > 0, x — > oo and y(x, 0) = 0, was given in [5]. The solution takes the form of the symmetrized 
(upper) curve in Fig. [H which represents an evolving grain boundary groove [24|. [5] . 

Alternatively, the above solution may be extended smoothly to the interior of an obtuse- 
angled wedge, as shown. These and the homothetic solutions in an acute-angled wedge, ob- 
tainable from the same type of reduction, make up the "open-angle" solutions [2D] . 
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Figure 1. Self-similar evolution by evaporation near a grain boundary groove (symmetrized upper curve) 
or equivalently by condensation in an obtuse-angled wedge. Output times are at t = 1 and t = 4. 



The homothetically shrinking simple closed-curve solution is the circle r = V2(to — r.) 1 / 2 . 
Self-intersecting closed-curve "flower- head" solutions were given in [T]. 

Steady-state solutions, invariant under time translations, are simply the straight lines. Like- 
wise, the solutions that are invariant under spatial translations in a particular direction are 
simply the straight lines in that direction. Another simple known solution is the travelling 
wave solution of the form y — ct = G(x), which is invariant under the translation in space and 
time generated by d/dt + cd/dy. This is the well-known "Calabi grim reaper" [17j . taking the 
appearance of a shepherd's crook 

y — ct = — log cos(cx). 
c 

Uniformly rotating solutions of the form r = G (eft— at) (r and (ft being plane polar coordinates) 
follow from reduction by the symmetry generated by a(—yd/dx + xd/dy) — d/dt. 

Expanding or contracting rotating solutions of the invariant form r/vt = G{(ft — ^ logt) or 
r/y/—i = G(cft — log(— t)) follow from reductions under the one-parameter group generated 
by a linear combination of generators of dilatations and rotations 

^^^^ddd^ d d 

r = r 4 + crr 5 , r 4 = x— + y— + 2t— , r 5 = -y— + x—. 

ox oy at ox oy 

A classification of the types of solutions of the associated reduced ordinary differential equa- 
tions, but without solving them explicitly and without referring to symmetries, was given re- 
cently by Halldorsson [18j . 

2.1 Reciprocal transformations 

It is well known that the class of nonlinear heat equations is stabilised by reciprocal transfor- 
mations which can sometimes be used to generate new solutions from old [22]. A reciprocal 
transformation may be viewed as a map of the graph of a solution u(x, t) of 



Ut = {D(u)u x ) x 
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to that of a solution u'(x',t') of 

u' t , = {D'(u')u' x ,) x , . 
The map (x, t, u) — > (x', t, v!) is defined by 

dx = udx + D(u)u x dt, dt' = dt, v! = — . (6) 

u 

It can be shown that this induces the transformation of diffusivities 

D'(s) = -U (-) . (7) 



s 2 \ s 



The diffusivity of interest here 
1 



D{u) 



l + u 2 

is seen to be invariant under the reciprocal transformation, 

D'{v!) = — ^ = D{u'). 

This raises the possibility of using the known solutions and generating new solutions of © and 
hence new curve flows by quadrature. However, it turns out that the reciprocal transformation 
of solutions of Q acts geometrically trivially on solutions of the curve shortening equation. 

Proposition 1. Suppose u(x,t) is a solution of ([3]) and y(x,t) the corresponding solution of 
the curve shortening equation. Let (x',t',u') be the reciprocal transformation of (x,t,u). Then 
the induced transformation on solutions of satisfies (x',t',y') = (y,t,x + c), c = const. That 
is, for each t, the reciprocal transformation of a known solution of (J3j) induces a reflection (up 
to additive constant) of the solution y(x,t) in the line y = x. 

Proof. Equation ([3]) and the reciprocal transformation ([6]) permits us to write 
x' = u(x\,t)dxi — / Vo(t\)dt, 

JO J a 

where a is a constant and Vo(i) = — D(u(0, t))u x (0, t). In terms of Cartesian coordinates (x,y), 
x' = y x { Xl ,t)dxi- Voih) = y(x,t)-y(0,t)+ / D(y x (Q,t x ))y xx (Q,t 1 )dt 1 

JO J a J a 

= y(x, t) - ?/(0, t)+ I y t {0, t x )dh = y(x, t) - y{0, a). 



Then 

u' = l/u = l/y x = x y = x x >. 

Now as y' satisfies y' x , = u' just as y x = u, then it follows that y' = x + Hit'). However, since 
y'(x',t') is to satisfy the curve shortening equation 

„,/ Vx'x' 



i + (?4) 2 ' 

it follows that H = const. 
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Let /i be a function such that the reflection of its graph defines a function fi- in general f\ 
and /2 will satisfy different differential equations and this can be very useful. For example, 
the nonlinear diffusion equation with diffusivity D(s) = s~ 2 , seen from ([7]) to be directly 
transformable to the classical linear heat equation (k = 1), is one of many integrable equations 
that can be linearized in its potential form by the hodograph transformation [12]: {x 1 , y') = (y, x) 
given in Proposition [TJ However, from a geometric point of view a curve and its reflection are 
indistiguishable. Thus invariance of © under reciprocal transformations does not increase the 
class of known curve flows. This highlights the difficulty of constructing solutions of the curve 
shortening equation and partly explains why so few interesting exact solutions are known. 



3 Functionally separable nonlinear heat equations 

Apart from symmetry methods, one can seek reductions that arise from second or higher order 
differential constraints rather than the first order differential constraints implied by classical Lie 
reduction. Unfortunately, there is no general procedure for seeking such reductions and much 
effort has gone into devising new reduction strategies for partial differential equations. Thus 
in Doyle and Vassiliou [14] the authors managed to classify all one-dimensional sourceless heat 
equations 

u t = {D(u)u x ) x (8) 

that admit separation of variables in some field variable u. That is, one asks for a change of 
field variable u = m(u), such that the image of ([8]) under the change of variable, constrained by 
the additively separable condition 

u xt = 0, 

is a differential system of finite type that can therefore be solved by ordinary differential equa- 
tions. Indeed the resulting differential system has the general form 

u t = f{u)u xx + g(u)u 2 x , u xt = 0, (9) 

for some functions / > 0, g. It is proven in [TJ] that for any such pair /, g there is a change of 
dependent variable that transforms ([9|)i to ([8]). It turns out that system ([9]) has a maximal 3- 
parameter solution space for any given pair /, g. In this manner the authors obtain exactly nine 
distinct diffusivities D up to the maximal transformation group that preserves the canonical 
form © for which the maximal 3-parameter solution space is achieved. In many cases the 3- 
parameter solution was constructed. This considerably extends the list of nonlinear diffusion 
equations for which explicit solutions are available. One equation on the Doyle-Vassiliou list is 
the nonlinear heat equation 

-(it?).- (10) 

Thus the differential 1-form 

u x 

uj = udx H 7T dt 

l + u A 

is closed on solutions of (|10p . It is easy to see that any function y(x,t) satisfying dy = uj is 
a solution of the curve shortening equation ([2]) . The solutions of (fTUI) constructed in [14] are the 
functions 



u(x, t) = U(a(x + a),a 2 (t + b)) , 
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where a, b and a^O are arbitrary constants and U is one of the functions 
U(x, t) = tan x, 



U(x,t) 
U(x,t) 



x 



V-x 2 -2t 
±1 

y/ e 2{x-t) _ I 



-rr, x sinhx 

U(x,t) = -== (11) 
V — cosh x — e a 



U(x,t) 

U(x,t) 
U{x,t) 



± cosh x 



V — sinh 2 x + e 24 

sinx 

V cos 2 x — e 2 * ' 

sinx 



V cos 2 x + e 



2i 



Each solution (jllh provides a 3-parameter family of explicit curve shortening flows except for 
the complex- valued, (fTT]) 1 , 

The function y(x,t) may be obtained from U(x,t) simply by integrating in x, then adding 
a suitable function of t. 

1- (|HP i then leads to the Calabi "grim reaper" travelling wave. 

2. (Illh o integrates to the well-known shrinking circle x 2 + y 2 = 2 (to — £)■ 

3. Integration of (IHP g merely produces a horizontal version of the vertical grim reaper. 
4- (|lip i is complex- valued, not considered further in the current practical context. 

5. (|TT]) £ is equivalent to (fTT|) 7 by a rotation in the xy-plane. 

6. The final two solutions are related to those previously presented by King [21] , expressed 
in the time-reversed form as examples of finger growth; these deserve closer inspection. 

3.1 Exact heat flow of a convex curve 

The studies of Gage-Hamilton [16] and Grayson [T7] on flow by curvature of embedded plane 
curves is a justly celebrated chapter in differential geometry. 

Theorem 1 (Gage-Hamilton). Let 70 : I x — > M 2 be a convex curve embedded in the plane. 
Let 70 evolve by the curve shortening flow. That is, 

— (x, t) = k(x, t)n(x, t), 7(x, 0) = 70 (x), 

where n(x,t) is the Euclidean curvature ofj(x,t). Then the curve remains convex and becomes 
circular as it shrinks in the sense that 

1) the ratio of the inscribed radius to the circumscribed radius approaches 1; 

2) the ratio of the maximum to the minimum curvature approaches 1 ; 

3) the higher order derivatives of the curvature converge to zero uniformly. 

(|ll|) e provides the only known example of an explicit, non-self similar curve shortening flow 
in case the initial curve 70 is a closed, convex embedded plane curve which is not a circle. We 
use results described in Section [21 The sixth function in (|11|) is 
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In the case of (|12p. the corresponding solution of ([2]) is 

y = t — In ^cos x + sj cos 2 x — e 2t ^j . (13) 

Clearly if y is a solution of ([2]) then so is —y. The two solutions join smoothly along y = and 
can be jointly expressed in the simple implicit form 

cosh y — e~ t cos x = 0. (14) 

This solution is recorded in [21] and shown here to arise from functional separation. Solutions of 
this type are also studied in [13] and referred to as 'Angenent ovals', although explicit solutions 
are not written down in the latter paper. We will now study some properties of this solution. For 
each t G (— oo,0) equation (114p defines a closed, convex "oval-shaped" curve which is symmetric 
about the x- and y-axes for 

x G (- cos" 1 (e*) , cos -1 (e*)) , y G (- cosh" 1 (e _i ) , cosh" 1 (e"*)) . 

By analogy with an ellipse, eccentricity may be defined as 



e(i) = \A ~ (^max/ymax) 2 



cos 1 (e*) 
cosh" 1 (e - *) 



2 



b\+o(\tf), 



showing approach to circularity (e — )■ 0) as t approaches the extinction time 0. The curvature 
at each point (x, t) of the curve is 



k{x, t) 



e 1 cos x 



Ve" 2 * - 1 



For each t G (— oo, 0), the maximum curvature occurs at x = with value K max = e */ V 'e 2t — 1, 
while the minimum occurs at the extremities along the minor axis, x = ± cos" 1 (e 4 ) with value 
Kmin = l/Ve" 2 * — 1. Hence 

lim max = lim e~* = 1, 

verifying the Gage-Hamilton theorem and more specifically showing that the ratio of maximum 
to minimum curvature is an exponential function converging to unity. Of course, the curvature 
itself is an unbounded function of time as the flow continues toward extinction. 

It transpires that for this curve flow, arclength can be expressed as a function of time in 
terms of the incomplete elliptic integral of the first kind. 

At very early times t < 0, the solution with x G (— 7r/2,7r/2) is asymptotic to two grim 
reapers joined smoothly and approaching each other with constant speed: 

y = ±(t- logcos(x) - log 2) + 0(e 4t sec 4 (x)). 

A question that may be posed is: why construct curve shortening flows by first solving (|10|) 
rather than solving ([2|) directly? Remarkably, it can be shown that 

Proposition 2. The image of the curve shortening equation ([2]) under the transformation y i— >■ 
m(y) does not have a (maximal) 3-parameter family of joint solutions with the linear wave 
equation y x t = in any field variable m(y) except m = identity in which case the solution gives 
rise to the grim-reaper flow. 
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-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

Figure 2. Evolution by heat shrinking flow of the curve defined by coshy — 5 cos x = 0. 

Proof. Construct the image of the curve shortening equation ([2]) under the change of variable 
y i — y m(y) and apply the Cartan-Kahler theorem to the differential system consisting of the 
transformed curve shortening equation and the constraint y x t = 0. ■ 

This is in sharp contrast to the rich separability properties of nonlinear heat equations (|8j) 
as discussed in Section El Thus starting with ([TO]) appears to be an important first step in 
constructing non-trivial curve flowsQ. 



3.2 Decaying periodic solution initially close to square wave 

For a curve fixed at two end-points, we prescribe the boundary conditions 

y = 0, X = 0,£. 

We now define dimensionless space and time variables X = x/£, Y = y/£, r = Bt/£ 2 , in terms 
of which the curve shortening equatiorH is 

Yr = (15) 

to be solved subject to boundary conditions 

y = 0, X = 0, 1 (16) 

and continuous initial conditions Y(X, 0) = Yq(X). 

Integrating the seventh member of the list (JTJJ and then applying translational and scaling 
invariance transformations, we obtain 



Y= K ln [ e X p(^[r-r o] ) ' ^ 



1 In view of the relationship between u(x,t) satisfying (| lOf) and y(x,t) satisfying j2), one might try for the 
higher order constraint y xxt = rather than y xt = 0. This is possible but the calculations rapidly become very 
complicated as the order increases. 

2 Also known in the materials science community as the Mullins equation [24] . 
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with K, Xq, to arbitrary constants. It may be verified that (|17p satisfies (|15p . In fact, when 
r — To is large and negative, 

In 2 1 / e -2K 2 [r -r] \ 

Y = — — K[t — t ] + — ]ncos(K[X - X ]) + O 



K K u,/ \cos(K[X - X Q ]) I ' 

This shows that away from the singularities of In cos {X — Xq), the solution is asymptotic in the 
distant past to the Calabi "grim reaper" solution of the curve-shortening flow |17j . At all times, 
the solution (|17p is a deformation of the grim reaper solution but now it is extended smoothly 
and periodically, without singularities over a domain of any length. Unlike in the grim reaper 
solution, there are fixed points so that we may apply the Dirichlet boundary conditions (|16p . 
which lead to 




exp(K^T-T J) / 



with 

The amplitude of Y(X, r), which is the value of |^(XT, t)| at X = ±n/(2K), is approximated by 

Y max = K[t - t] + Qn2)/K + O (exp (-K 2 \t - To|)) . 
In terms of dimensional quantities, 

y max = ^^"^ + i— + O (exp (-nVSfe, " t]/i 2 )) , (19) 

where y max = ^max and to = £ 2 t~o/B. This shows that at early times, the amplitude decreases 
linearly as a function of time. 

As an example, the solution is graphed for the case n = 5 and with the minus sign preceding 
the right hand side of (|18p . The solution at early times is shown in Figs. [31 HI Although for 
large to, the initial condition resembles a periodic square wave, it actually converges pointwise 
to a differentiable bounded and periodic grim reaper as to approaches oo. 

When K 2 [t — to] is large and positive, the solution is smoother, approximated by a single 
sine wave, 

Y = ± e M-K 2 [r-T }) s . m{Rx) + Q (exp( _ 2if 2 [T _ To])) ) (2Q) 

whose amplitude decreases exponentially in time. This is the sinusoidal solution of the classical 
linear diffusion equation that approximates the curvature driven diffusion equation at large 
times in the small-slope approximation. 

The number of extrema n within the fixed domain [0, 1] may be freely chosen. For large 
values of the parameter -KT 2 to, as in Fig. 02 the initial condition is described well as a periodic 
square wave, resembling a diffraction grating. The time scale for decay may be viewed as the 
time at which the formal expression (I19p for Y max is zero. This is 

£ 2 ln2 
* — *o H — 2 2p - 

Fig. [3] evidences considerable smoothing of the initial conditions but at this time, the solution 
does not yet resemble a simple sinusoid. At larger times, the solution is close to a single sinusoidal 
wave. In the regime of the sinusoidal profile, from (|20p decay times are shorter, of the order 
of the standard time (£/nir) 2 /B that is familiar from exponential decay of Fourier modes in 
linear diffusion models (e.g. [TU]) when the half-wavelength tjn is the typical distance between 
neighboring regions of high and low mass concentration. 
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Figure 3. Filled curve: exact solution with t — to = —0.07, K = 5n. Dashed lines: approximate bounds 
y = ±K[T -T] + (ln2)/K. 
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Figure 4. Filled curve: exact solution with t — tq = —0.02, K = 5ir. Dashed lines: approximate bounds 
y = ±K[T -T] + (hi2)/K. 

4 Surface evolution by evaporation and condensation 



For some metals such as gold, surface diffusion persists for several thousand years as the dom- 
inant surface transport mechanism but for others such as magnesium, after a time less than 
a day, surface evolution occurs predominantly by evaporation-condensation. As described by 
Mullins [23], the Gibbs-Thomson relation for evaporation rate in terms of deficit from equi- 
librium pressure over a curved surface, leads to a second-order equation for diffusion by mean 
curvature. In terms of two-dimensional Cartesian coordinates (x, y) and time t, the Mullins 
equation for points on a material surface is 



Vt = B 



Dxx 
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where B is constant. This equation applies to two dimensional cross sections of solids when 
surface nano-scale features such as grooves, ridges and furrows extend rectilinearly into the 
third dimension. Because of the nonlinearity, very few useful exact solutions to this equation 
are known [8], except in a linear approximation. Two decades ago [5], the exact solution was 
constructed in parametric integral form for nonlinear surface evolution near a symmetric grain 
boundary, with constant slope at the grain boundary groove, initial flatness and zero displace- 
ment at infinity. A similar procedure produces the more general "open angle" solutions for 
deposition in a wedge [20] , Subsequently, the fourth-order Mullins equation for curvature-driven 
surface diffusion on an almost-isotropic material, was solved with boundary conditions repre- 
senting a symmetric grain boundary [6l [28] [?]. 

In the context of surface evaporation, the parallel asymptotes of the grim reaper solution 
represent a long thin metallic foil that is evaporating at the ends. The constant travelling-wave 
speed of the grim reaper solution shows that evaporation will take place at a constant rate 
that depends on foil thickness as well as the evaporation coefficient B. The thickness Ax is the 
distance between the two asymptotes, Ax = ttB/c. Hence the steady evaporation rate at the end 
of a strip of metallic foil will be c = ttB/Ax. For example, for a foil of a few microns in thickness 
made of the volatile metal Mg, this rate will be of the order of one millimetre per millenium. 

In the context of metal surface smoothing, the periodic solution of the previous section 
predicts the smoothing of initial conditions that resemble a diffraction grating. For a surface 
fixed at two end-points, we prescribe the boundary conditions 

^ = 0, x = 0,£. 

Physically, this corresponds to the surface being clamped and shielded from the surrounding 
atmosphere outside of the exposed spatial domain [Q,£]. 

4.1 Evaporation from anisotropic crystals 

C. Herring [19], showed that for anisotropic crystals, surface energy is proportional to b(4>)('y + 
j"(4>)), where <j) is the polar angle arctany x , 7 is surface tension and b is an evaporation 
coefficient. An equation of the form 

Vt = D(y x )y xx , (21) 

implying the standard nonlinear diffusion equation ([8]) with u = y x , may be regarded as an 
anisotropic form of the Mullins evaporation-condensation equation 



with surface slope-dependent anisotropy factor B{y x ) = D{y x )(l + y x ), which originates in the 
physical derivation from a constant multiple of b{(p){^ + 7" (</>)) (e.g. |27]). Crystalline materials 
are indeed anisotropic, with the evaporation coefficient minimized when the cut surface is aligned 
with crystal planes. The group of equivalence transformations of this class of equations includes 
the general linear group GL(2,M). By the polar decomposition theorem, each invertible linear 
transformation can be decomposed as an othogonal transformation followed by multiplication by 
a positive definite symmetric dilatation matrix. Under a rotation about the origin by angle a, 

x = x cos q — y sin a, y = x sin a + y cos a, i = i, 




u + tan a 



sin a + y x cos a 



Vxx 



u = y x = 



1 — u tan a 



Vx = 



cos a + y x sin a 
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The axis x = has been rotated by angle —a/2. Writing 6 = arctany s , 
n • ey _ cos 6 _ 1 1 

yt = — vi = — 75 — \K = — -i i — i — Vt- 

ne y cos{0 — a) cos a 1 + y 5 tan a 
Hence, by rotation, the isotropic nonlinear diffusion equation (|21|) is equivalent to 

= D (yx)yxx, D(s) = -. — : ■ 



(s sin a + cos a) 2 \ s sin a + cos a 

Tritscher [26] used this device of rotational equivalence classes to solve integrable forms of 
fourth-order surface diffusion equations. In particular, after rotation by angle a = vr/2, 

1 ( -\ 

D(m) = ^d{^- 

Vs V Vx 

By following the tt/2 rotation by a trivial reflection x —> —x, we recover the result of the 
reciprocal transformation ((7|). 

By the principal axis theorem, a positive symmetric matrix can be written as O t QO, where 
O is a rotation matrix and Q is diagonal, = with aj > 0. Therefore, to consider the 
effect of an additional dilatation, we need only consider the effect of a diagonal rescaling: 

2 ( a i - \ 

x = aix, y = a 2 y, y t = a-JJ — y x y xx , 

which although trivial, allows us to construct simple anisotropic models from the rotationally 
invariant isotropic model. 



5 Anisotropic analogs of isotropic model solutions 

The scale invariance group still applies to the general anisotropic diffusion equation (|2ip . There- 
fore both expanding and shrinking types of similarity solution exist but they may be significantly 
different from those of the isotropic model. The grain boundary groove solution still exists, as 
can be seen from the solvability of the reduced boundary value problem on M + x M + , 

u = F( P ), P = xr 1 /\ 

^F>(p) = ±[D(F)F>(p)], 

F(0) = m > 0, F(p) ->0, p^oo. 

The grain-boundary groove solutions for all models have some common features. Although it 
represents a highly anisotropic material, the linear model groove solution approximates that of 
an isotropic model for groove slopes of up to 0.5, which was used in Mullins' original paper [24] , 
However, whereas the linear model predicts that the groove depth increases in proportion to m, 
the groove depth increases very slowly, of order (logm) ' 5 at large m for the isotropic model [3]. 



5.1 Anisotropic homothetically shrinking closed curve 

For anisotropic evaporating materials, the closed-curve homothetic solution represents the fixed- 
shape cross section of an evaporating wire, which is circular when the material is isotropic. 
The homothetic evaporating closed-curve solution satisfies 



y = [t - t]V 2 G(p), p = x[t -t]- 1 / 2 , 
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— + H -G\p) = G"{p)D{G'), 

G"(0) = 0, G(po) = 0, G"(p)^oo, p^po- 
This implies 

u = F(p), F = G', 

P -F'{p)-D'{F)(F'f = F"{p)D(F), (22) 

F(0) = 0, F(p) -> oo, p -> po- 

For some functions D(u), the homothetic closed-curve solution does not exist. For example, 
with the linear model with constant D, the general solution satisfying F(0) = must be 

F(p) = ^Aerf{p) = A [ P e~ s2/i ds, AeR. 
2 Jo 

Since this cannot take an infinite value at any point p = po G R, the closed-curve homothetic 
solution does not exist. 

Let us now assume that F(p) ~ Aq(pq — p) v . If a homothetic closed-curve solution exists, 
G(po) G (—oo,0), and 

rpo 

= G(p ) = G(0) + / F{s)ds. 

Jo 

This can hold only if the integral does not diverge, implying 
-1< v < 0. 

Now we suppose that at large u,D(u) is asymptotic to a power law D{u) ~ Dqu u . Then by 
balancing terms at the leading order in po — p, (|22|) implies 

(i) n = 1/V, which is less than —1, and po = — 2vDqA^ u , or, 
(ii) n = which is less than —2. 

For the linear model, the anisotropy factor is B(u) = (1 + u 2 )D which diverges when the 
curve is vertical. In the current application, we are interested in cases of realistic anisotropy for 
which the evaporation coefficient and the anisotropy factor are bounded, and the latter with 
a minimum value greater than zero: 

VueM, < B < B{u) < < oo. 

That statement must be true, independent of the orientation of the coordinate axes. For example 
after rotation by ir/2 it must be true that D(u) — > Dq G (0, oo) as u — > 0. By rotating back to 
the original orientation, this implies 

D(u) = u- 2 D(u) ~ Dou' 2 . 

This physical restriction rules out scenario (ii) in the above leading-order analysis. The simplest 
example satisfying the restriction is D(u) = 1/(1 + (flu) 2 ) (/3 G R constant) for which the 
anisotropy factor B(y x ) varies between 1 and f3 2 . In this case the homothetic solution is simply 
an ellipse 

? 2 

2[t -t] + 2f3 2 [t -t] = 
elongated in the direction of weakest evaporation. 
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5.2 Grim reaper solution for anisotropic material 

The grim reaper solution is a travelling wave constrained between two vertical asymptotes. This 
implies a steady state solution of the equation ut = [D(u)u x ] x . For arbitrary D(u), there exists 
a two-parameter steady state solution 

u s (x) = K _1 (cx + c 2 ), 

where 

K{u) = [ D(s)ds 
Jo 

which is an increasing invertible function because D(s) > 0. This integrates formally to 

y(x,t) = y (t) +/ u s {xi,t)dx 1 . 
Jo 

It then follows from (|21|) that yo(t) can only be a constant-velocity translation, yo = c(t — to). 
This travelling wave solution does not necessarily have a vertical asymptote. For example, there 
is no such asymptote when D is constant. However if D(s) satisfies the physical requirement 
-D(s) ~ s~ 2 , it follows that K(u) has a finite limit as u — > oo, therefore K -1 has a vertical 
asymptote at some location x = xq, duplicated at x = —xq if D(s) is an even function, as is 
commonly the case when the cc-axis denotes the orientation of the crystal planes from where 
evaporation is weakest. 



6 Anisotropic models allowing functional separation 

The classification of Doyle and Vassiliou gives all functions D(s) for which functional separation 
of variables is possible in the general form 

u = f(u) = v(x) + w(t), (23) 

with / an invertible function. 

The simplest example of a physically feasible isotropic model is simply that obtained from 
the isotropic diffusion equation by unequally rescaling x and y. For example, Fig. [2] could be 
dilated in one direction, displaying a non-homothetic closed curve approaching a homothetically 
shrinking ellipse. 

The only other member of the Doyle-Vassiliou list with realistic anisotropy is much more 
complicated: 

D(u) = Dq cos(z(^4u)), Au = / (cos s y 3/2 ds, -it/2 < z < it/2. (24) 

J o 

This model is close to isotropic when A = \[2. In that case, it is easy to show that when 
u (= Vx) is small, B{u) = 1 + 0{u A ) and B — > 2 as u — > oo. In Fig. [U the function D(u)/Dq for 
this weakly anisotropic model is compared to D{u)/D$ of the isotropic model. 

For the sake of completeness, we construct the special solution u(x,t) compatible with (123|) . 
that was not given explicitly by Doyle and Vassiliou [13]. The parameter A may be changed by 
rescaling u. For convenience, without loss of generality we now set A to 1. Also we may set Do 
to 1 by using Dot as the time coordinate. From the general approach of Doyle and Vassiliou |14j . 
u is a sum of separated functions v{x) and w(t) satisfying 



v'"{x) + -[v' (*)] 3 = 0, 



(25) 
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Figure 5. Nonlinear diffusivity function D{u) for isotropic evaporation (lower) and anisotropic model 
(upper). 



w'(t) = W(w) = cosw 



v cos v + - (v ) sin v 



+ sinu> 



-v smv + — (v ) cosv 



By construction, the following are first integrals of (|25p . which can be verified by substitution: 
A = v"(x) cosv + -^{v') 2 sinv, fj, = —v"(x) sinv + T^iv') 2 cost). 

It then follows that w'(t) = A cosw + (isinw, and that 
2 f dw 1 f v ds 



I 



R 2 J sm(w + 5) 



R Jo y/cos(s - 5) 



(26) 



with constants xq, to, R = ±2 1 / 2 (/u 2 + A 2 ) 1 / 4 and 8 = arctan(A///), whose arbitrariness is of 
little consequence. Let v'(x) = W(v), so v"{x) = WW'(v). Then 



WW'{v) cos v + -W 2 sinv = A, 
1 nsin 2 u 1, 



-WW'{v) sin v + - VF 2 cos v = (jl 



-Atanu + -VF + -W z cosv = u 

2 cos v 2 

VF 2 = 2(/x cos f + A sin t>) = R 2 cos(v — 5) 

1 f v ds 



'/(x) = ±R^ cos(v — 5), 



R Jo ycos^s - 5) 



When we include the parameter Bo, the general solution for v and w in terms of elementary 
functions , the standard elliptic integral F{9 | |) and the standard Jacobi elliptic function sn (X | | ) , 
is 



2 arctan e R Bot / 2 — S and v = 2 arcsin f —= sn ( — - — ^ ° 

.V2 V y/2 



+ 5, 



where 



V2. 

-R 1 V 



-^-F (arcsin ( V2sm 
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X 

Figure 6. Portion of solution curve for anisotropic model with R = 1 and 5 = 0, at t = 0. 



3.5 - 

3 - 



2.5 - 




Figure 7. Portion of solution curve for anisotropic model at t = —2.0. 

Since u is a function of w(t) + v(x), we see that the choice of parameters R and 5 has little 
consequence on the form of the solution u(x,t). A shift in the phase variable 5 has the effect 
of a translation in x by xq(5,R). A change of amplitude R has the effect of rescaling x and t 
to Rx and i? 2 £. 

Modulo a time-dependent vertical translation, the curve y(x, t) is obtained from u(x, t) by 
integration. Since the integrands in (fMj) and (|26l) must be real valued, the constructed solution 
has a truncated domain. An example is given in Figs. [6H3 

As t approaches — oo, the solution approaches a steady state. Since D(u) is symmetric, this 
steady state is a symmetric grim reaper. It is approximated in Fig. [8] by taking t = —8.0. 

The curve is not symmetric, as can be seen in Fig. [71 This asymmetry causes its domain to 
shift slightly from right to left. For t < 0, the curve has a vertical asymptote at moving location 

rTT/2-w(t) 

x = x u = I (cos s) ' ds 

Jo 
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Figure 8. Portion of solution curve for anisotropic model at t = —8.0. 



and a local minimum at 



X = X£ = — 



Mi) 

I (cos s) 






7 Conclusion 

Of the exact solutions to the curve shortening equation known to us, most can be obtained by 
Lie point symmetry reductions. The two interesting solutions that cannot be constructed in this 
way, can indeed be recovered by functional separation of variables for the standard nonlinear 
diffusion equation ([3]) that is obtained from the curve shortening equation by differentiation. 
The classification obtained by Doyle and Vassiliou [H] of nonlinear diffusion equations that 
admit functional separation of variables, leads to the two exact non-self-similar solutions with 
non-trivial initial conditions that appear to be achievable in an approximate sense in applica- 
tions. In addition, it leads to a new separated solution for a physically realistic anisotropic 
evaporation-condensation diffusion equation. Although the second order nonlinear surface evo- 
lution equations for slope u(x,t) (= y x (x,t) admit a number of possibilities for functional se- 
paration of variables in Cartesian coordinates, we have proved that this is not possible for the 
equation ([2]) governing y(x,t), nor is it possible in a coordinate system consisting of canonical 
variables for a symmetry other than translation. This is in contrast with the point symmetry 
analysis, which leads to a richer array of possibilities for the evolution of y(x,t) than for the 
evolution of u(x,t). 

The invariance of the isotropic equation (I10p under the well known reciprocal transformation 
was shown (Proposition [1]) to lead to no new planar curve heat flows. The group of geometric 
equivalence transformations of the class of general anisotropic equations (121j) includes not only 
the reciprocal transformation in the guise of a reflection in the plane, but the whole general linear 
group. The equivalence is shown by carrying out the equivalence transformations explicitly. 
Physical evaporation coefficients must have a positive real value when the surface is oriented 
along the crystal planes. Since physical restrictions must be independent of orientation of the 
coordinate axes, it follows from the equivalence transformations that the nonlinear diffusivity 
D(y x ) must behave like y~ 2 at large-slope. This also allows for the existence of a closed non- 
circular homothetic solution which cannot exist unless D{u) decreases faster than it . The 
Doyle- Vassiliou classification produces another anisotropic model that satisfies this physical 
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requirement. An exact solution has been constructed, involving Jacobi elliptic functions and 
other inverse integrals of rational functions. 

Exact solutions sometimes have the advantage of leading to concise conceptually simple 
relationships. For example, Fig. [3] demonstrates the efficacy of a simple expression for wave 
amplitude of a corrugated nano-scale surface in the early stages of smoothing by evaporation- 
condensation when the system cannot be adequately described by a linear model. However, 
exact solutions can be obtained only in very special cases of initial and boundary conditions, so 
that approximate numerical solution methods will continue to be important. 

Acknowledgements 

This paper is submitted in appreciation of the valuable on-going contributions of Professor 
Willard Miller Jr. The first author gratefully acknowledges support by the Australian Research 
Council under project DP1095044. 

References 

[1] Abresch U., Langer J. The normalized curve shortening flow and homothetic solutions, J. Differential Geom. 
23 (1986), 175-196. 

[2] Angenent S., On the formation of singularities in the curve shortening flow, J. Differential Geom. 33 (1991), 
601-633. 

[3] Arrigo D.J., Broadbridge P., Tritscher P., Karciga Y., The depth of a steep evaporating grain boundary 
groove: application of comparison theorems, Math. Comput. Modelling 25 (1997), no. 10, 1-8. 

[4] Bluman G.W., Reid G.J., Kumei S., New classes of symmetries for partial differential equations, 
| J. Math. Phys] 29 (1988), 806-811, Erratum, fT Math. Phys.\ 29 (1988), 2320. 

[5] Broadbridge P., Exact solvability of the Mullins nonlinear diffusion model of groove development, 
| J. Math. Phys] 30 (1989), 1648-1651. 

[6] Broadbridge P., Tritscher P., An integrable fourth-order nonlinear evolution equation applied to thermal 
grooving of metal surfaces, \TMA J. Appl. Math.\ 53 (1994), 249-265. 

[7] Broadbridge P., Goard J.M., Temperature-dependent surface diffusion near a grain boundary, 
| J. Engrg. Math] 66 (2010), 87-102. 

[8] Cahn J.W., Taylor J.E., Overview no. 113 surface motion by surface diffusion, Acta Metall. Mater. 42 
(1994), 1045-1063. 

[9] Cao F., Geometric curve evolution and image processing, Lecture Notes in Mathematics, Vol. 1805, Springer- 
Verlag, Berlin, 2003. 

[10] Carslaw H.S., Jaeger J.C. Conduction of heat in solids, 2nd. ed., Clarendon Press, Oxford, 1959. 

[11] Chou K.-S., Li G.-X., Optimal systems and invariant solutions for the curve shortening problem, Comm. 
Anal. Geom. 10 (2002), 241-274. 

[12] Clarkson P.A., Fokas A.S., Ablowitz M.J., Hodograph transformations of linearizable partial differential 
equations, \SIAM J. Appl. Math.\ 49 (1989), 1188-1209. 

[13] Daskalopoulos P., Hamilton R., Sesum N., Classification of compact ancient solutions to the curve shortening 
flow, J. Differential Geom. 84 (2010), 455-465. la7Xiv:0806. 17571 

[14] Doyle P.W., Vassiliou P.J., Separation of variables for the 1-dimensional non-linear diffusion equation, 
Unternat. J. Non-Linear Mech. 1 33 (1998), 315-326. 

[15] Galaktionov V.A., Dorodnitsyn V.A., Elenin G.G., Kurdyumov S.P., Samarskii A. A., A quasilinear equation 
of heat conduction with a source: peaking, localization, symmetry, exact solutions, asymptotic behavior, 
structures, J. Soviet Math. 41 (1988), 1222-1292. 

[16] Gage M., Hamilton R., The heat equation shrinking convex plane curves, J. Differential Geom. 23 (1986), 
69-96. 

[17] Grayson M. 1987 The heat equation shrinks embedded plane curves to round points, J. Differential Geom. 
26 (1987), 285-314. 



Symmetry and Separation in Surface and Curve Flows 



19 



[18] Halldorsson H.P., Self-similar solutions to the curve shortening flow, arXiv: 1007. 1617 

[19] Herring C, Surface tension as a motivation for sintering, in The Physics of Powder Metallurgy, Editor 
W.E. Kingston, McGraw-Hill, 1951, 143-179. 

[20] Ishimura N., Curvature evolution of plane curves with prescribed opening angle, Bull. Austral. Math. Soc. 
52 (1995), 287-296. 

[21] King J. R., Emerging areas of mathematical modelling, ! ii. Soc. Lond. Philos. Trans. Ser. A M ath. Phy s. Eng. Sci.\ 
358 (2000), no. 1765, 3-19. 

[22] Kingston J.G., Rogers C, Reciprocal Backlund transformations of conservation laws, \Phys. Lett. A\ 92 
(1982), 261-264. 

[23] Malladi R., Sethian J. A., Image processing via level set curvature flow, Proc. Nat. Acad. Sci. USA 92 (1995), 
7046-7050. 

[24] Mullins W.W., Theory of thermal grooving, fj. Appl. Phys] 28 (1957), 333-339. 

[25] Olver P.J., Sapiro G., Tannenbaum A., Invariant geometric evolutions of surfaces and volumetric smoothing, 
\SIAM J. Appl. Math.\ 57 (1997), 176-194. 

[26] Tritscher P., Integrable nonlinear evolution equations applied to solidification and surface redistribution, 
PhD Thesis, University of Wollongong, 1995. 

[27] Tritscher P., An integrable fourth-order nonlinear evolution equation applied to surface redistribution due 
to capillarity. 17. Austral. Math. Soc. Ser. 51 38 (1997), 518-541. 

[28] Tritscher P., Broadbridge P., Grain boundary grooving by surface diffusion: an analytic nonlinear model 
for a symmetric groove, \Proc. Roy. Soc. Lond. A\ 450 (1995), no. 1940, 569-587. 



